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Abstract. A useful sampling-reconstruction model should be sta- 
ble with respect to different kind of small perturbations, regardless 
whether they result from jitter, measurement errors, or simply from 
a small change in the model assumptions. In this paper we prove 
this result for a large class of sampling models. We define differ- 
ent classes of perturbations and quantify the robustness of a model 
with respect to them. We also use the theory of localized frames to 
study the frame algorithm for recovering the original signal from 
its samples. 

1. Introduction 

The sampling and reconstruction problem includes devising efficient 
methods for representing a signal (function) in terms of a discrete (finite 
or countable) set of its samples (values) and reconstructing the original 
signal from the samples (see e.g., PQ El El El EZl E2] and the reference 
therein). In this paper we consider a very general sampling model where 
the signal is assumed to belong to a finitely generated shift invariant 
space and the sampling is performed on an irregular separated set and 
is averaged by finite Borel measures. The main focus of this paper is on 
describing and quantifying "admissible" perturbations of the sampling 
model which may result from altering the sampling set (jitter) (see e.g. 
P [71 Ej), or the averaging sampling measures (measuring devices) or 
the generators of the underlying shift-invariant space (see e.g., [51 [TS]). 

As recently became customary in sampling theory (see e.g. [U [21 
[TTl [T9l [20l [2n [22]). we mesh operator theory techniques and those of 
shift invariant and Wiener amalgam spaces |[13j. The latter provide 
us with relatively straight-forward proofs while the former allow us to 
keep in sight our objective. In section 2 we show that all the properties 
of our sampling model can be encoded in the sampling operator U . 
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The sampling model admits reconstruction if its sampling operator is 
bounded both above and below. Our first goal is to show that any 
and all of the small perturbations mentioned above result in a small 
perturbation of U in the operator norm. This will prove the stability 
of sampling in our model with respect to those perturbations and the 
corresponding estimates we obtain will quantify this stability. Our 
second goal is to show how a frame algorithm can be used to reconstruct 
signals in our sampling model. Finally, our last goal is to show that the 
reconstruction error due to the perturbations we describe is controlled 
continuously by the perturbation errors. 

The paper is organized as follows. In section 2 we describe our 
sampling model, introduce relevant notions and notation, and cite a 
few preliminary results. The main results are presented in section 3. 
Perturbation results addressing our first goal are in subsection 3.1. 
There we prove that a set of sampling remains such under a small 
perturbation of the sampling measures and/or the generators of the 
shift invariant space. It is also shown that sampling remains stable with 
respect to a perturbation of the sampling set itself. In subsection 3.2 
we show that, in case of a signal in a Hilbert space, a frame algorithm 
can be used to reconstruct the function from its samples. We also use 
the results of the previous subsection and the theory of localized frames 
to show that under mild additional assumptions a set of sampling for 
a Hilbert shift invariant space is also a set of sampling for a chain of 
Banach shift invariant spaces to which the frame algorithm extends. 
In subsection 3.3 we study the dependence of the reconstruction error 
upon the perturbation errors. The proofs of the results in section 3 are 
relegated to section 4. 

2. Description of the sampling model 

This section is primarily devoted to introduction of the sampling 
model we use in this paper. We also present most of the necessary 
notation and cite some of the preliminary results that will be used 
later. 

The signals we are studying in this paper are represented by functions 
/ G L^(M'^), for some p E [1, oo] and d eN. Moreover, we assume that 
/ belongs to a shift invariant space 

(2.1) VP{<^) = {Y, Cl^k : C e (F(Z'^))M}. 

Here $ = (0^, . . . , (jfY ^ vector of functions, $fc = $(• — k), and C — 
(c^, . . . , d'Y is a vector of sequences belonging to {(P{'L'^))^'^\ Among 
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the equivalent norms in (£P(Z'^))(^') we choose 

r 

11^*11 (^p(Zd)){r) = ^ ||c*||£P(Zd). 

i=l 

In order to avoid convergence issues in (12.11) we assume that the set 
{(f)^{- — k), . . . , (f)^{- — k);k G Z"^} generates an unconditional basis for 
VP{^). In particular, we require that there exist constants < nip < 
Mp < oo, such that 

(2.2) 

mp\\C\\^enz^)),r> < II E Cl'^kh. < Mp||C||(,.(^.))„, VC G {F{Z') f\ 

The unconditional basis assumption (12.21) implies [3] that the space 
VP{(^) is a closed subspace of U'{R'^). 

Since we are interested in sampling in V^{^) we add an assump- 
tion that would make all the functions in these spaces continuous and, 
therefore, pointwise evaluations will be meaningful. To this end, we as- 
sume that all generators $ belong to a Wiener-amalgam space ( Wq)^^^ 
as defined below. For 1 < p < oo, a measurable function / belongs to 
WP if it satisfies 

(2.3) ||/|k.= f5^esssup|/(a; + A;)|A < oo. 

If p = oo, a measurable function / belongs to 1^°° if it satisfies 

(2.4) 11/11 woo = sup{esssup \ f{x + k)\} < oo. 

Hence, W°° coincides with L°°(R'^). It is well known that are 
Banach spaces [13], and clearly Wp C Lp. By (iyp)('-) we denote the 
space of vectors ^ = ('0^, • • • , ^^)^ of VK^-f unctions with the norm 

r 

ll^ll(W'P)('-) = ll"^! WP- 

i=l 

The closed subspace of (vectors of) continuous functions in Wp (re- 
spectively, (WpY"-^) will be denoted by (or {W^Y'^). 

In this paper we are interested in average sampling performed by 
a vector of measures. We denote by M{R'^) = Mo(R'^) the Banach 
space of finite complex Borel measures on R''. The norm on Ai{M.'^) is 
given by \\fj,\\ = J^^ d\fi\{y), i.e., the total variation of a measure By 
{A4{W^)Y^^ we denote the space of vectors 'Ji = (yU^, . . . , yU*) of measures 
from AdlW^) with the norm ||7^||(>i(]Rd))(t) = X]j=i \\f^''\\- The symbols 
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Ms{R'^) i{MsiR'^)y*'>), < s < oo, will be used for the subspace of 
mIm.'^) ((A^(M^))W) of all (vectors of) measures /x G M(M.'^) such that 
(1 + |a;|)" e L\R'^,d\fx\), i.e., /(I + \x\yd\fi\{x) < oo. By ModM'^) 
{{Moo{R'^)Y^^) we denote the space of all (vectors of) measures with 
compact support. Clearly Al5(]R'^) C Mr(R^) for < r < s < cxd. 

For fi G Ai{R'^) and a measurable function on M!^, the convolution 
of the function (p and the measure /i is defined by 

(0*/i)(x)= / (P{x-y)dfi{y), xGM'^. 

When we have a vector of measurable functions $ = (0^, . . . , 0^)^ and 
a vector of finite complex Borel measures 7? = (z^^, • • • then the 
convolution $ * 7^ is the r x t matrix given by 

(0^ * /i^ ... 0^ * /i* 
: : 
(f)^ * fl^ ... 0'" * /i* 

Let J be a countable index set and X = {xj : j G J} be a subset 
of M°'. The reconstruction problem in our sampling model consists of 
finding the function / G V'^($) from the knowledge of its samples 

(/ * jt)iX) = {(/ * Jt)ix,) = {if * /.i)(x,), ...,(/* /x*)(a;,))},ej. 

When t = 1 and /i = 5o, i.e., /i is the Dirac measure on R'' concen- 
trated at zero, then {f*~jl){X) = {fixj)}j^j and we obtain the classical 
(ideal) sampling model. When d'Jx = "^dx, where \E' G (L^(R'^))*^*) and 
dx is the Lebesgue measure on M'^, i.e., 'Jx is absolutely continuous with 
respect to the Lebesgue measure, then we write (/ * '^){X) instead of 
(/ * ^iid our model is reduced to the case analyzed in [5]. 

Definition 2.1. Let 1 < p < oo and X = {xj : j G J} be a countable 
subset of W^. We say that X is a set of sampling for V^{^) and 'Jx (or, 
simply, a /i -sampling set for V^{^)) if there exist constants < Ap < 
Bp < oo such that 

(2.5) ApWfU. < ||(/*7?)WII(£.(j))W < BpWfU., for all / G V^^). 

If d'Jx = "^dx then a 7?-sampling set X will be called a \E'-sampling 
set and, if t = 1 and fi = Sq, then X will be called an ideal sampling 
set. To ensure that an upper bound Bp in (12.51) always exists (see (14. 2p ) 
we restrict our attention only to separated sets X. 

Definition 2.2. We say that X is separated if there exists 6 > such 
that inf ij^j^ijLj \xi — Xj\ > 6. The number 5 is called the separation 
constant of the set X. 
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It is not hard to extend our results to the case of a finite union of 
separated sets. We do not, however, pursue this relatively trivial but 
space consuming generalization. 

Definition 2.3. Let 7? G {M{R'^)Y*\ $ E {W^lY'''^ satisfy (E2]), and 
X = {xj,j G J} C M"^ be a separated set. The sampling model is the 
triple {X, <l>, 7?)- The sampling model (X, $, /?) is called p-stable if X 
is a 7^-sampling set for p G [1, oo]. 

Given a sampling model (X, $, /?) we proceed to define its samphng 
operator. 

Definition 2.4. The sampling operator U = U(^x,<i>,~il) '■ (i^ (1^^))^'^^ — > 
(£P(J))(*) is defined by UC = {f * 7? )(X), where / = E C^^k G 

We can think of U t X r matrix of operators 

/ f/i'i . . . W'^ 

£/= ; ; 

V f/i'* . . . f/'"'* 

where for each 1 < i < r and I < I < t the operator f/*'' is defined by 
an infinite matrix with entries {U^''')j,k = {4>^*fJ'^){xj — k), j E J, k E Z"^. 
The operator norm of U is given by ||t/||p^op = Sz=i X]I=i 

The following proposition shows that all the interesting properties 
of a sampling model (X, $, 7?) are, indeed, encoded in the sampling 
operator U. The proof of this result follows immediately from (12.21) 
and fl23D . 

Proposition 2.1. T/ie sampling model (X, $, 7?) is p-stable if and only 
if there exist < r]p < Pp < oo such that for all C G the 
sampling operator U satisfies 

(2.6) '7p||C'||(^p(z'*))('-) < II^C'||(^p(j))(t) < /3p||C||(£p(2d))(r). 

The next lemma is, essentially, a nutshell for many of the results in 
this paper. 

Lemma 2.2. Let (X, $, 7?) be a p-stable sampling model and U be its 
sampling operator satisfying l{2.6\) . Let also (X, 9, 'a) be a sampling 
model such that its sampling operator satisfies \\U — Uj\\\ < rjp. 
Then (X, B, 'a) is also p-stable. 

Proof Let C G {i^ {Z'^))'^'-\ Then 

\\UaC\\ + \\uc\\ 

< — f^All ||C'||(£p(z<*))('-) + Pp\\C\\(ep(Z''))(^)- 
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Therefore, since \\U — U^\\ < rjp, then we have 

(2.7) \\UaC\\ 
On the other hand, since 

'7pl|C'll(£p{Zd))M < l|f^C'll{<?p(J)){*) < - f^A)C'||(^p(j))(t) + ||f/AC||(^p(j)){t) 

Hence, 

(2.8) (r^, -\\U-Ua\\) < ||f/AC||(,,(^))W. 

Since — t/A|| < %, the conclusion of the lemma follows from (12. 7p . 
(ESD, and Proposition O □ 

3. Main Results 
In this section we collect the main results of our paper. 

3.1. Admissible perturbations of a sampling model. 

In practice, shift invariant spaces are used to model classes of signals 
that can occur (or that are allowed) in applications. However often, 
the functions in a shift invariant space model only give approximations 
to the signals of interest. For this reason, we begin with a result where 
the perturbation of a sampling model is due to a small change of the 
genetators of the underlying shift invariant space. 

Theorem 3.1. Let (X, $, 7?) be a p-stable sampling model for some 
p G [l,oo]. Then there exists eo > such that the sampling model 
(X, 0, 7?) is also p-stable, whenever Q G {Wq)^'^'> and ||$ — 9||(;yi)(r) < 
eo- 

The above result means that if 7? e {M{W^))^'\ $ G ( Wo^)(") satisfies 
(12.21) . X = {xj,j G J} C M"^ is a separated -sampling set for \^p($), 
and G satisfies the assumptions of the theorem, then there exist < 
< Bp < CO such that 

(3.1) A;||(?|Up<||((7*7?)(X)||(,p(^))(.)<i?;i|^7|Up, forall^7G\/P(e). 

In the proof of this result in section 4 we will provide explicit estimates 
for eo and the bounds and Bp. 

As a consequence of Theorem 13.11 we have the following results that 
were first proved in [5]. The proofs now are immediate: we apply The- 
orem 13.11 with dji = "^dx for Corollary 13.21 and 'Jl = So for Corollary 
1331 
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Corollary 3.2. Let ^ e {L\W^)Y^\ $ e (Wf^Y''^ satisfy HJj, and 
X = {xj,j E J} C M.'^ be a separated "i/ -sampling set for V^{^). Then 
there exists eo > such that X is a "^-sampling set for ^^(6), whenever 
e G (^^o')^"^ (^rid ||$-e||(vi.i)w < eo. 

Corollary 3.3. Let $ G {WqY''^ satisfying ^ES) and X = {xj,j G 
J} C 6e a separated ideal set of sampling for Then there ex- 

ists eo > such that X is an ideal set of sampling for V^{Q), whenever 
e G {WflY'^ and 11$ -e||(^i,i)M < e < eo. 

In practice, signal samples are obtained using measuring devices with 
characteristics that are not fully known, and the measurements reflect 
local averages rather than exact sample values. Thus, a sampling mea- 
sure 7^ is a model that approximate the characterisitics of a measuring 
device. For this reason, the next theorem describes the case when the 
perturbation is due to some uncertainty about the characteristics of 
the measuring devices, that is a perturbation of the vector of measures 
7?. 

Theorem 3.4. Let (X, $, 7?) be a p-stable sampling model for some 
p G [l,C)o]. Then there exists eg > such that the sampling model 
(X, $, "a') is also p-stable, whenever!^ G {Ai(W^)Y''^^ and 

Again, if X, 'Jx, 'a, and $ satisfy the assumptions of the theorem 
then there exist < < Bp < oo such that 

(3.2) A'jfU. < IK/* < 5;||/|U., for all / G V^{<l>), 

and the explicit estimates for eo, A'p and will be given in section 4. 

Considering and 'a in Theorem 13.41 such that dji = ^dx and 
d'a = Tdx we obtain the following direct corollary (see also [5], Theorem 
3.3]). 

Corollary 3.5. Let ^ G (^^(R'^))^, $ G {W^Y''^ satisfy and 
X be a separated "^-sampling set for V^^($). Then there exists eo > 
such that X is a T -sampling set for V^{^), whenever V G {L^{W^)Y*^ 
and 11^ - r||(^i(jjd)){t) < eo. 

As a consequence of Theorems 13.11 and 13.41 we obtain the follow- 
ing combined perturbation result and its corollary, which is essentially 
Theorem 3.4 in [5]. 

Theorem 3.6. Let (X, $, 7?) be a p-stable sampling model for some 
p G [l,oo]. Then there exists eo > such that the sampling model 
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{X,Q,'a) IS also p-stable, whenever 'a G {M{R'^)Y^\ 9 G {W^^Y^'^ 
and 11$ — 6||(iyi)(r) + \\ jj, — a ||(yK(Rd))(t) < ^o- 

Corollary 3.7. Let ^ G {L\R'^)Y^\ $ G {W^^Y''^ satisfy HJj, and 
X be a separated "if -sampling set for Then there exists eo > 

such that X is a V-sampling set for V^{Q), whenever V G {L^{W^)Y*\ 
e G {W^Y''^ and ||<l> - e||(iyi)M + 11^ - r||(ii(igd))(t) < eo. 

An error in the location of the sampling points {xj} is what is often 
called jitter error (see e.g., [6l[7j and the references therein). This error 
can be modeled as a perturbation of the sampling set X. For this 
reason, our next perturbation results deal with an altered sampling set 
X = X + A = {xj + 5j}ieJ, where A = C W^. We use the 

standard notation for ||A||^ = sup{||5j|| : j G J}. 

Theorem 3.8. Let (X, $, 7?) be a p-stable sampling model for some 
p G [l,C)o]. Then there exists eo > such that the sampling model 
{X + A, $, 7?) is also p-stable, whenever || A||^ < eo- 

Remark 3.1. The above theorem is an analog of Theorem 3.6 in [6J, 
where r = t = 1, p = 2, and fi = fi^ = Sq. 

As a direct corollary of Theorems 13.61 and 13.81 we get the following 
combined result. 

Theorem 3.9. Let (X, $, 7?) be a p-stable sampling model for some 
p G [l,C)o]. Then there exists eo > such that the sampling model 
(X+A,e,^) zs also p-stable, whenever^ G {M{R'^)Y^\ 6 G {W^Y''^ 
^"^d ||A||^ + 11$ - 6||(n/i){r) + 117? - "« ||(A4{Md))W < eo- 

We leave it to the reader to formulate other perturbation theorems 
resulting from different combinations of Theorems 13.11 [37il and l3.8[ We 
conclude this section with a slightly stronger version (due to Lemma 
of Theorem EH 

Theorem 3.10. Let (X, $, 7?) be a p-stable sampling model for some 
p G [l,oo] and U be its sampling operator. Let also (X + A, 9, "c?) 
be a perturbed sampling model with the sampling operator Ua- Then 
for every e > there exists eo > such that \\U — Ua\\ < e, whenever 
'a G (7W(R'^))W, e G (1^0^)^"^ and 

ll^lloo + ll*^ " + 117? - ll(A^(]Rd))(t) < eo. 

3.2. Perfect reconstruction and localized frames. 

In this section we show that a frame algorithm can be used to recon- 
struct / G V^^($) from its samples. We also obtain a useful modification 
of the above results using the theory of localized frames developed in 
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[T5] (see Definition [33]). In the previous section, the number p E [l,oo] 
was fixed, that is, we stated, for example, that if X is a 7?-sampling 
set for V"^($), then X is a /i -sampling set for V^{<d) for the same 
p G [1,00], as soon as 6 is sufficiently close to $ in the appropriate 
norm. Here, we claim that if X is a 7^-sampling set for V^^($), then X 
is a 7^-sampling set for V^i^Q) for all p G [1, 00], as soon as 9 is suffi- 
ciently close to $, $ satisfies a mild decay condition, and Jt belongs to 
7Vls(I^'^) for some s > d. It is natural to ask whether one can replace 
in the above statement with for some q G [1, 00]. Under 

certain assumptions the answer is "yes" , but it turns out to be a much 
harder problem as shown in [2]. 

Definition 3.1. Let be a Hilbert space of functions and V a closed 
subspace of H. Let {"^Xj = ii^lj^ ■ ■ ■ ^4'ij)'^}jeJ be a countable collec- 
tion of vectors of functions in V. We say that {\I'^^}jgj is a frame for 
V if there exist constants < A < B < 00 such that 

AWfWn < < BWfWn, for all feV, 

where (/, ^x,) = ((/, , • • • , (/,<)) e C*. 

Remark 3.2. Notice that the above is not quite the standard definition 
of a frame in a Hilbert space. This is due to the way we defined 
the norm in ([2]). Nevertheless, it is easily seen that {"^x^yjej is a 
frame for V according to the above definition if and only if {"0* ., i = 
1,2, ... ,t, j G J} is a frame for V according to the standard definition. 
The frame bounds, however, may be different. 

Definition 3.2. Let \^ be a closed subspace of the Hilbert space Ti. 
Let {"^xj = (^Ij j • • • y'4'ij)^}jeJ be a frame for V. The frame operator 
associated with the frame {\E'^^}jgj is the operator S : V —>■ V defined 
by S{f) = Eiej(/'^^.)*^.' all f e V. The (canonical) dual 
frame {"^x of the frame {"^x is a sequence of vectors given by 
= i^l^,. . where rx, = ^'V^,, l<s<t. 

Remark 3.3. It is well know that a frame operator 5* is bounded, invert- 
ible, self-adjoint, and positive [12]. Hence, the canonical dual frame is 
well defined. There may exist other dual frames but we will refrain 
from defining the notion. 

The next proposition shows that a frame algorithm can be used to 
reconstruct a function from its samples. 

Proposition 3.11. Let $ G ( Wo^)('^^ 7? G {M{R'^)Y^\ and X be a 7?- 
sampling set for y^($). Then there exists a sequence of vectors of func- 
tions {'^Xj}jej, which is a frame for V^{^) and (/, \I'^.^) = (/ * 'jt){xj) 
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for all f G V'^{^) and j G J. Moreover, every function f G can 
be recovered from the sequence of its samples {(/ * 'Ji){xj)}j<^j via 

(3.3) /(a;) = ^(/*7?)(a:,)^.,(x), 

where is the dual frame of o^nd the series liS. 3\) con- 

verges unconditionally in y^($). 

The frame {^^.^ jjgj constructed in the previous proposition will be 
called a {ji ^ X)- sampling frame for V^^($). The main idea of this sec- 
tion is to use the fact that if such a frame is localized then it is also a 
Banach frame [TB] for \^^($), p G [1, oo). 

Remark 3.4. Observe that, in general, the frame operator S is the 
product of the analysis operator T : V — * {P{J))^^\ defined by Tf = 
{(/, ^x,)}ieJ = {((/, , •••,(/, V'x,) )}jGJ and its adjoint, that is ^ = 
T*T. Since $ generates a Riesz basis, it is immediate that in case of 
a (7?, -^)- sampling frame its analysis operator is isomorphic to the 
sampling operator U = f/(x.$,7?)- 

Definition 3.3. Let V" be a closed subspace of the Hilbert space 
H. Let {"^Xj = ii'lj, ■ ■ ■ be a frame for V, and {Gk = 

{g\, ■ ■ ■ , gk)'^}k&d be a Riesz basis for V, i.e., a condition similar to 
(12.21) is satisfied. We say that the frame is (polynomially) 

s-localized with respect to the Riesz basis {Gk}k&<i) if 

(3.4) KG,,C,)|<Ci(l + |x,-A;|)-, 
and 

(3.5) |(G,,C,)l<C2(l + |a:,-A:|)-^ 

for all j & J and k G Z"'. Here, the constants Ci,C2 > are inde- 
pendent of J and k, |(Gfc,^y| = EI=i ELi I (^L I ' {^fcjfcez'' is 
the dual Riesz basis of {Gk}ke:Z'i, and \ {Gk, is defined similarly to 
\{Gk,%)\. 

Remark 3.5. Let K be a closed subspace of a Hilbert space Ti. As- 
sume that {Gk = {gl, ■ ■ ■ ,gkY}k&d is a Riesz basis for V. The dual 
Riesz basis of the Riesz basis {Gk}kezd is the sequence of vectors 
{Gk = (gl, . . . ,gl)'^}k(zzd satisfying {Gk,Gf) = 6kil, where I is the 
r X r identity matrix, and 6ki is the Kronecker delta. Since a Riesz 
basis {Gk} is also a frame, {Gk} is, in fact, the canonical dual frame 
for {Gk}- In this case it is the unique dual frame. 
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Definition 3.4. Let $ = . . . , 0'^)^ G (^^o^)^'^) C (L2(M^))M and 
s > d. We say that $ is an s-localized Riesz generator for \^^($), 
denoted $ G W^, if 

• {$fc = — fc)}fceZ'' generates a Riesz basis for \^^($), i.e., 
condition fl2.2p holds for p = 2; 

• The components of $ satisfy the decay condition 

(3.6) \<P\x)\ < C',{1 + \x\r , 

for all 1 < i < r and some Cq > independent of x G M.'^. 

Remark 3.6. If $ G Ws, then (12. 2p holds for every p G [1, oo] as shown 
in [3]. 

The following is the main result of subsection 3.2. 

Theorem 3.12. Let s > d, <^ e Ws, and Jt G (A1^(M'^))W. A ssume 
that X is a 'Ji -sampling set for V'^{^), and is the {ji^X)- 

sampling frame for V"^($) . Then 

• X is a 'Ji -sampling set for V^^($) for all p G [1, oo]. 

• If {"^xj} is the dual frame for xj}jeJ J then 

(3.7) / = ^(/ * 7?)(x,)$.^, for all f G V^($), 

where the series converges unconditionally in V"^($), p G [1, oo). 

Next, we combine Theorem 13.121 with the perturbation results of the 
previous section. The proofs are immediate. 

Theorem 3.13. Let s > d, <i> e Ws, and 7? G {MsiM.'^)Y^\ A ssume 
that X is a separated 'Jl -sampling set for y^($). Then there exists 
eo > such that for every 9 G Wg satisfying 11*^* — 0||(v^i)('-) < eo; there 
exists a (jt , X) -sampling frame {\E'^^.}jgj for V^{<d). Moreover, 

• X is a 'Jx -sampling set for ^^(0) for all p G [1, oo]. 

• If {"^Xj} is the dual frame for {^Xj}j<^j, then 

f = J2^f * 7?)(^,)^.,, for all f G V^iQ), 

where the series converges unconditionally in V^{Q), p G [1, oo). 

Theorem 3.14. Let s > d, ^ E Ws, and jt G (A<s(K°'))W. .4s- 
sume that X is a separated 'Jt -sampling set for V"^($). Then there 
exists eo > such that for every Ix G (A1s(R"'))*^*^ satisfying ||7? — 
"c?||(_yvi(Rd))W < ^0; there exists an (a , X) -sampling frame {'i!x^}ji^j for 
V^^($). Moreover, 
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• X is an Tx -sampling set for V''{^) for all p G [l,C)o]. 

• If {"ifxj} is the dual frame for then 

/ = * for all f e 

where the series converges unconditionally in V^{^), p G [1, oo). 

Theorem 3.15. Let s > d, <^ e Ws, and 7? G {MsiM.'^)Y^\ A ssume 
that X is a separated Ji -sampling set for Then there exists 

eo > such that for every A = {Sj, j G J} satisfying ||A||^ < eo there 
exists a {j2 ^ X + A)-sampling frame {^xj}jeJ fof V^i^). Moreover, 

• X + A is a ~j2 -sampling set for for all p G [1, oo]. 

• If{'i>xj} is the dual frame for xj}jeJ) then 

f = J2if * 7?)(^. + Sj)^^,^ for all f G V^i^), 

where the series converges unconditionally in p G [1, oo). 

Theorem 3.16. Let s > d, ^ e Ws, and 7? G {Ms(M'^)Y^\ Assume 
that X is a separated Ji -sampling set for Then there exists 

Co > such that for every Q G Ws and~a G (7Vfs(R'^))'^*'' satisfying ||$— 
0||(VKi)('-) + IIT'* ~ "^II(.A/i(]R'*))(*) < ^0; there exists an {'a , X) -sampling 
frame {'^xj}jeJ for V^{Q). Moreover, 

• X is an 'a -sampling set for ^^(0) for all p G [1, oo]. 

• If {"^xj} is the dual frame for {^^^j^gj, then 

f = J2if * « )(^.)*.,, for all f G y^(e), 

where the series converges unconditionally in V^{Q), p G [1, oo). 

Theorem 3.17. Let s > d, ^ e Ws, and 'Jl G {Ms{^'^)Y^^ ■ Assume 
that X is a separated 'Ji -sampling set for V^($) . Then there exists eo > 
such that for every A = {6j, j e J}, Q e Ws, and "a G {Ms{M.'^)Y^^ 
satisfying ||A||^ + ||$-©||(^^i)(r) + ||7t'-"a ||(^(]Kci))(t) < eo, there exists 
an (a, X + A)-sampling frame {^xj}jeJ for ¥"^(0). Moreover, 

• X + A is an 7x -sampling set for ^^(0) for all p G [1, oo]. 

• //{^-r^.} is the dual frame for {^xj}jeJ, then 

f = Y^if * -S){x^ + 5,)^,., for all f G V^{Q), 
where the series converges unconditionally in Vp{Q), p G [1, oo). 
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Remark 3.7. The crucial result for the proof of the theorems in this 
section is Jaffard's non-commutative extension of the classical Wiener's 
Tauberian Lemma (see Theorem 5 in [I5]). It states that if an invert- 
ible matrix has an off-diagonal decay defined by inequalities similar to 
(13.41) and (13.51) . then the inverse matrix has the same off-diagonal de- 
cay. There exist other extensions of Wiener's Lemma which deal with 
different types of off-diagonal decay (see, for example, [10l[T6]). Many 
of those could be used to obtain results similar to, say. Theorem 13.171 

3.3. Imperfect reconstruction. 

In practice, we know that a perturbation exists because of imperfec- 
tions of measuring devices, errors, etc. However, we can only estimate 
this perturbation and may not even know its nature. Here we show 
that even if we use a model (X, $, 7?) for reconstructing a signal from 
a perturbed model (X, G, 7?) (or vice versa), the reconstruction error 
depends continuously on the perturbation in the cases studied above. 

As before, let U be the sampling operator for a p-stable sampling 
model (X, and Ua be the sampling operator for a perturbed 

model (X, 6, "c?), where X = X + A = {xj + The sampling 

operator Ua can be thought of as a t x r matrix of operators given by 



Ua 



\U'^ 



U'a' ) 



where for each 1 < z < r and 1 < / < t the operator f/^' is defined by 
a bi-infinite matrix with entries (f/^')j,fc = (^* * c^^[xj + 5j — k), j G J, 
keZ'^. 

We let U* be an operator defined by the following r x t matrix of 
operators from (£p(J))(*) into (£p(Z^))('): 



/ f/i.i . . . f/i.* \ 



U* 



\ . . . W'^ J 



where for each 1 < z < r and 1 < / < t, the operator is defined by 
a bi-infinite matrix with entries (t/*'')j,fe = (0* * ~ where z 



denotes the conjugate of the complex number z. The operator (Ua)* is 
defined similarly. Notice that this definition implies ||f/*||p^op = ||t/||p,op, 
and (U*)* = U. Moreover, if U satisfies (12. 6p . then U* satisfies 



(3.^ 



Vp\\D\\(ep{j))W ^ ||f^*-D||(fP(zd))('-) < Pp\\D\\{ep(j))it): 
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for all D e {P'{J)Y^\ Observe also that if p = 2 then U* is, indeed, 
the Hilbert adjoint of U . Hence, if the sampling model (X, $, 7?) is 
2-stable, U*U is isomorphic to the frame operator 5* for the sampling 
frame {\E'a;^.}, see Remark |3.4[ Therefore, U*U is invertible and positive. 
Moreover, the operator {U*U)~^U is a left inverse for the sampling 
operator U and it is isomorphic to the synthesis operator used for the 
reconstruction. Hence, the importance of the following result. 

Theorem 3.18. Let (X, $, 7?) be a p-stable sampling model for some 
p G [1, oo]. Assume that its sampling operator U satisfies Ii2.6\) and the 
operator U*U is invertible. Let e G (0, —(3p + ^/ Pp + Vp) ^'^^ 
be a perturbed sampling model such that its sampling operator sat- 
isfies \\U — Ua\\ < e- Define v = u^e) = r]~'^e{e + 2[3p) . ThenO <1, 
the operator U^Ua is invertible, and 



\\{U*U)-'U* - iUlUA)-'Ul\\ < ^ e + 



Vp \ ^-^ 



Remark 3.8. Observe that if p = 2 we do not need to require invert- 
ibility of U*U . As we mentioned above, it follows automatically. 

Remark 3.9. If in Theorem 13 . 1 8 1 we let r = t = 1, p = 2, and /x = /i^ = 
5o, then we obtain an analog of Theorem 3.3 in [B]. 

Let (X, $, 7? ) be a p-stable sampling model for some p G [l,C)o]. 
Assume that its sampling operator U satisfies (12. 6p and the operator 
U*U is invertible. We define the reconstruction operator R = -R(x,<i),7?) • 
(^P(J))W ^ VP{^) by 

RD=J2m*U)-'U*D]m.-k), 

D = {d\...,d'Y in {^P{J))^\ 

Then as an immediate consequence of Theorems 13.101 and I3.18^ we 
have the following result. 

Theorem 3.19. Let (X, $, 7?) be a p-stable sampling model for some 
p G [l,C)o]. Assume that its sampling operator U is such that U*U is 
invertible. Let R be the reconstruction operator. Then for every e > 
there exists eo > such that for every A = {Sj, j G J} , Q E (WqY'^^ 
and 'a G (A^(]R'^))*^*^ satisfying 

ll^lloo + ll*^ ^ 0|l(VFi)M + 117? - "^ll(A^(Rd))(*) < eo, 

we have 

\\R{{g * 7?)(X + A)) -f\\^^<e,f=J2 Cl<^k, 9 = CjO^, 
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for die e 

Theorem 13.191 tells us that the reconstruction error is, indeed, con- 
trolled in a continuous fashion by each and all of the perturbation errors 
studied in this paper. 

Our final result is a combination of the above theorem with the 
results of section 3.2. 

Theorem 3.20. Let (X, $, 7?) be a 2-stable sampling model such that 
$ G Ws and 'Ji E (^^^(IR'*))''*^ Let R be the reconstruction operator 
for (X, $, 7?). Then for every e > there exists eo > such that for 
every A = {6j, j G J}, G Ws, and "a G satisfying 

ll^lloo + 11*^ ^ 0|l(VFi)M + 117? ~ "^ll(A^(Rd))(*) < Co, 

we have 

\mg * 7?)(X + A)) -f\\^^<e,f=J2 Cl<^k, 9 = ^ C^Qk, 

for allpe [l,oo] and all C e {e'{Z'^) f'\ 

The proofs in the following section show implicitly how numerical 
estimates for eo in Theorems 13.191 and 13.201 may be obtained. 

4. Proofs 

4.1. Auxiliary results. 

We begin with technical results that are needed for the main proofs. 

Lemma 4.1. Let (p e Wq , and /i G M(M'^). Then: 

(4.1) (f)* fi e Wq, and 

(4.2) <2'^||0||^^i||^||. 

Proof. Note that if yU = 0, the proof is immediate. Assume now yU 7^ 0, 
i.e. ||/i|| > 0. Let e > be given. Since G Wq, then is uniformly 
continuous in M'^. Therefore, there exists 6 = 6{e) > such that 

(4.3) \4>{w) — (j){wi)\ < - — -, whenever — < 6. 
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Let zq G M*^ be given, and let 2; e M"^ be such that — ^oll < ^- Then 
we have 

\{^* fx){z) - {(p* fx){zo)\ = / (j){z - y)dfi{y) - (f){zo - y)dfi{y) 



< 



{(p{z -y)- (pizo - y))d^{y) 
\(p{z-y) - (l){zQ-y)\d\ii\{y). 



Since || {z — y) — {zq — y)\\ = \\z — zq\\ < 6, for all y G M'^, then it follows 
from gis]) that J^^ \{(p{z -y) - (j){zo -y))\d\fi\{y) < J^, ^^d\fi\{y) = e. 
Since zo and e > are arbitrary, we obtain the continuity of * in 

Let us show (g^D- Let G and 12 e M{R'^) be given. Then 



■ * ^Wiyi = esssup 



J (j){x + k-y) dfi{y) 



< 



J esssup \(p{x + k - y)\d\ij\{y) < 

/ I J2 esssup + k-y)\ \ d\fx\{y) = J ||0(- - y)\\ d\ii\{y). 
.k&d ze[o,i]'* / 



Since ||0(- < 2^11011 vFi, for all y G we get 

w^d\^i\{y)< ! 2''U\\w-d\^l\{y) = 2^ 
Therefore, we get (14. 2p . 



□ 



The next proposition collects basic facts about Wiener amalgam 
spaces, shift invariant spaces V^{^), and separated sets in M"^. 

Proposition 4.2. Let $ G {Wf^Y'\ 7? G (;V/1(M'^))W, / = E C^^k, 

w/iere C G (^^(Z^))^, an(i <l>fe = $(■ - A;), for all k G Z"^. Let also 
X = {xj,j & J} be a separated set in M*^ with a separation constant 
6>0. Then 



(4.4) 
(4.5) 
(4.6) 
(4.7) 



$*7? G (l^o')^'''*^ 

1"^ * 7^ll(lVl)('->^') ^ 2'^||$||(jyl)(r) ||7^||(X(Rd))(t); 

yp($) c 1^0". /o^ a// 1 < p < 00; 



< \\c\\ 
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(4.8) \\f{X)U^,j)<U\\f\\w.. where U = U{5,p,d) = + lY'^. 

Proof. First, Lemma [4.11 immediately implies (14. 4p . 

Next, to prove (14. 5p consider $ = (0\ . . . , (j/Y' ^ ( W^o^)'" and 7? = 
(/iS . . . , /i*) G (;V/1(M'^))W. Then using we obtain 

t r 

II* * 7^ll(H/i)('-xt) = 5]^^ 110' */^^||iyi < 

i=l i=l 

t T 

Y,Y.'^'\W'MM\ = 2l$IUi)«||7^||(>i(M.))(*). 
j=l i=l 

Next, we prove (14. 7p . Consider 1 < p < oo and f = Yl ^I^k- For 

each 1 < s < r let a'^{l) = esssup |0'^(x + /)|, for all I G Z*^. Then 

a'e[o,i]'' 

Ikll^Hz;'') = ll'/''IUi- Consequently, ||a||(^i(2d))M = ||$||(vi/i)W, where 
a = (a^, . . . , a^)"^, and $ = (0^, . . . , 0'')"^. Hence, 

r r 

esssup I /(x + /) I <X15Z |c'(fc)|esssup|0^(a; + /-A;)| = ^(a^* |c^|)(/). 

^£[0,1]'* ,=1 ^^^d ^e[o,i]'* ,=1 

By using Young and triangular inequalities, we have 

r r 

ll/lkp < 5Z ||a^ * |c1l|£p < 5^ ||a1|£i||c^||£p. 

S = l 3 = 1 

Consequently, ||/||iyp < ||C||(£p(2d))r ||<l>||(^yi)(r). 

Next, let us show (gl]). Let / G Vp{^) be given. Then / = 
EfcezdC'J<l>fc, for some C G (£p(Z"'))(^'). Since (HTD implies / G Wp, 
it remains to show the continuity of /. Let us first consider the case 
1 < p < oo. We observe that C = L°°(M^) (see Theorem 2.1 
in [1]), and, hence, 

(4.9) ||/||Loo(R.)<rfl||/|kp, 

for some di > independent of /. Let /„ = X]|A;|<n ^I^k be a partial 
sum of /. Since $ G ( Wq^)*-'^\ then {/„}„eN is a sequence of continuous 
functions, and from (14.71) and (14. 9 p we obtain 



11/ ~ /n||L°°(R'*) < '^l||'^ll(VFi)('-) 
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Therefore, the sequence of continuous functions {fn}nm converges uni- 
formly to the function /. Thus, / is a continuous function as well. To 
treat the case p = oo, we choose a sequence {$n}n>i of continuous func- 
tions with compact support (see Theorem 3.1 in [1] for details) such 
that ||<l>„ - <l>||(iyi)(r-) ^ as n ^ oo. Set fn{x) = Ylik&'i ^I^nix - k). 
Since the sum is locally finite, then each is continuous. By using 
(14.71) once again, we estimate 

ll/n " /IIl°°(IR'*) ^ '^l||C'||(^oo(2d))(r) — $||(J^l)(r). 

It follows that the sequence of continuous functions {/„}n>i converges 
uniformly to /. Hence, / is a continuous function as well. 

Finally, let us prove (14.81) . Since X = {xj,j G J} C M is sep- 
arated with a separation constant 6 > 0, then infj^^ |xj — Xk\ > S. 
Consequently, there exist at most ([^] + l)'' sampling points in every 
(i-dimensional hypercube [0, l]'^ + Z, Z G Z''. Therefore, 

l/(a:,)r < (riv^+l)^esssup|/(x)r, 
and, hence, ||/(X) < ^\\f\\wp, for all / G W^, where Af = 

Using (14.41) and (14. 5p . we obtain the following result. 

Corollary 4.3. Let A : (Wo^)^^) x (A<(R'^))W — > (VF(i)(''^*) be defined 
by A($, 7?) = $*7^. Then A is a bounded bilinear form, and \\A\\ < 2'^, 
where 

\\A\\ = sup{||A(<l>,7?)||(^^^i)(.xt) : ||$||(vv^i)M < 1, ||7?||(A^(Rd))W < 1}. 

The following lemma, proved, for example, in [5J, states that a small 
perturbation of a Riesz basic sequence remains a Riesz basic sequence. 

Lemma 4.4. Let $ G (M^^)^'') satisfy Ii2.2\) . Then there exists eo > 
such that every 6 G (W^Y"^^ satisfying ||$ — B||(jyi)(r) < e < eo, also 
satisfies i\2.S\) . for some < m'p < < oo and 

(4.10) mp>mp — e and < \\^\\(^^n^(r) + e. 

4.2. Proofs for Section 3.1. 

Now we are ready to prove the first of our main results. 

Proof of Theorem 13.11 

Proof Assume that Ji G {M{R'^)Y^\ $ G {Wf^Y''^ satisfies (Q, and 
X = {xj,j G J} C M'^ satisfies f l231) . We want to find eo > such 
that whenever ||$ — 0||(WAi)(r) < e < eo, then (13.11) takes place for some 
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< Ap < Bp < oo. Assume < e < mp. Then, by Lemma WM 9 G 
{W^Y''^ satisfies (E^D and we can use representations g = J^kez^ ^I^k 
and / = Y.k&z^ Cl^k, C e )('■). Consequently, we have 



< He'll (£p(zd)){r) < — II y2, ^k^kWip — — 



< 



< 



+ 



LP 



k&d 



nip 
rrir, 



ll(/*?)WII(. 



\k&<l 



(£P(J))(*) 



A, 



-1 t 



rrir. 



1=1 



\k&<i 



IP (J) 



A-l * 



m 



1=1 



\k&<i 



iP{J) 



'-\\{9*1^) Wll (^p(j))(*) ' 



where 
(4.11) 



t. 



Since $ and 6 are elements of {W^Y''^ and 7? G (A<(M"'))(*\ then by 
n . we have S' = ($ - 6) * /x' G (^yo^)(''\ for Z = 1, . . . , t. Hence, 
using fl4.5p . (14.61) and condition (12.21) for g = J^kez^Cj^Qk, we have 

^1=1 \\i^k&'i ^k ^k)i-^)\\iP(^j) ^ 

2''A/'||C||(^p(2d))(r) 11$ - 6||(tyl)(r) ||y||(A^(Md))(i) < 



Therefore, 



M' 
p 



LP 



< 



+ 



Ap ^2'^Af\\^ - e\\ II 7^11 (M(R'i))(^) 



mprrip 



* l^)iX)\\{ep{j))M- 



LP 



+ 
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Hence, 

'ApTUp 2^M\\^ - 0||(l4^i)M||7^||(A^(Rd)){t) 

LP 



(4.12) V K K 



< 



||(^* 7?)(x)||(^p(j)){t). 



On the other hand, since 9 G ( Wo^)^'^^ and Jt G {M{M'^))^^\ it follows 
from g2D that {9^ * fi^ , . . . ,9"- * fi^) e { W^Y'"\ l = l,...,t. Therefore, 
f l47l) . dMD and the first of the estimates in (HTTOjl imply that 

< Ar||((J]cJe,)*7^)||(i^.)„ 

< 2'^A/'|| 7?||(^(iRd)){t)|| '^ClQk\\(wp)(-) 



yd 



< 



l|C|| l|0|Ui)« 

2'^-^ll7ril(_^(Kd)){t) 



< \IU\\LP- 



Hence, 
(4.13) 

/i)(x)||(^p(j))w < ' _ ^ — I hWip- 

\ Ulp t 

Using the estimates (14.101) and the left hand side of the inequality 
(I4.12p . we can obtain an explicit upper bound eo for e from 

ApTUp 2^A/'||7i ||(^(igd))(f) ^ _ ^ 



e nip — e 



This is equivalent to the quadratic equation 

e^ + C e- ^P""' 

" 2'^Ar||7/||(^(Kd))(t) 



where 



ApTUp 
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Let eo be the positive solution of the previous equation, i.e., 



1 / / 2 4:Apml 
Then, for < e < cq < rrip, we use ( KT^ . (gUSD, and fHTTU]) to obtain 

^/ _ ApUlp _ 2^A/'||7i ||(;^(iRd))(t) ^ 

^ ||*||(VFi)M + e nip-e 

, 2'^A/'||7?||(>!(Rd)){t)(||$||(w'i)M +e) 

= ) 

mp — e 

and the proof is complete. □ 
Proof of Theorem 13.41 

Proof. Let / = Cj$fc G C G (£p(Z'^) )('■). We have 

< IK/ * (7? - 7?)) (X) II + IK/ * 7?) (X) II 

= E IK/ * - + IK/ * « ) Wll(..(j)w 

= E ll((E cT*^^) * - WIK^(^) + IK/ * «)(^)ll(.^(.))(')- 

Since /i and a are in (A^(M'^))W, and $ G {W^^Y^'K then Proposition 
lO implies fi' = (0^ * (/i' - a^) , . . . , (p'' * (/i' - a')) G (W^o^)^'"^ for 
/ = 1, . . . , t. Using Proposition 14.21 once again we have: 

t 

ApWfWLP < X]A/'||C||(^p(2d))M||f]'||(^^i)(r) + IK/ * "a)(X)||(^p(j))(t) 
1=1 

< 2'^A/'||C||(^p(2d))M||$||(iyi)(r)||7? - "« ||(A^(R<i))(t) + IK/ * ■a)(X)||(^p(j))(t). 

Taking into account $ G (W^)^'^^ also satisfies fl2.2p . and / satisfies 
fl23D . then it follows 

„d»n|^|| ii^ii II-. , 2^Ar||/|Up||$||(H/i)M||7?-^ll(MR<*))W 



rup 



Hence, 
(4.14) 

. 2'^Ar||'^ll(W^i)w||7?-^ll(A^(K''))W 1 II / ||//-,--Vx.N|| 

\\j\\lp< |K/*a)(A)||(^p(j))(t). 



rUp 
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On the other hand, since / G V"^(<l>) satisfies (12.51) . we have 

IK/ * < IK/ * (7? - 7?))(x)||(,,(^))w + IK/ * 7?)(^)ll(..(j))w 



t 



= E IK/ * - + IK/ * 7?)(^)ll(..(j))w 

< E iK(E ^^'^'^^ * - wii^n^) + ^pIi/iu^ 



< 2''A/'||C||(£p(2d))(r)||$||(tyl)(r)||7? - "a ||(_^(Kd))(t) +i?p||/||LP. 

Using condition (12. 2p . we obtain: 
(4.15) 

^ , / 2^Ar||<|)||(^^i)M||7? - "a ||(^(i;d))W D V|.|| 
ll(/*tt)WII(£P(j))W < ^ — +Bp II/IUp. 



rup 



From (I4.14P and (I4.15p . by choosing 



ApUip 

eo - 



2'^Ar||$||(^^i)„' 
we obtain for < e < cq, 

" nip 

2'^AA||$||(^i)(.) 
-°« = -Dp H f ■ 



□ 



Proof of theorem 13.61 

The conclusion of the theorem is essentially obvious at this point. 
We proceed with a formal proof in order to obtain estimates for eo and 
the bounds and of X as an "c?-sampling set for ^^(0). 

Proof. Let < < i {^j^l^^^^^, ' c)) , where 



^ ll(A4{R'^))(*) 



Then, by Theorem 13. ![ X is a yU -sampling set for V^{Q) as soon as 
Moreover, 

^pII^IIlp <\\{g* 7?)(^)ll(^p(j))w < KWahr', for all g G v^^(e). 
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where 

_ ApUip 2'^7V|| /i ||(_A4(Rd))(t) 

and 

„ 2'^J\f\\ n ||(A^(]Rd))(t)(||$||(w'i)W + ei) 

B = . 

mp-€i 

Assume now that 

v4p(mp-ei) 

< 62 < 



2''Ar(||$||(j^i)M +ei)' 
Then, by Theorem 13 .4^ X is an "c^-samphng set for V^{Q) as soon as 

ll*^ - 0|l(tyi)M < ei and - "a ||(^(Kd))(t) < £2. 
Hence, if < e < cq = min{ei, £2}, we obtain the samphng bounds 
_ , 2'^AA(||$||(t^i)M+60 
mp - ei 

and 

2'^Ar(||$||(^,)M+ei) 

^p = -Dp ^ ^2, 

as soon as 

11$ - 0||(iyi)M + 117? - ll(A^(iRd)){t) < e < eo- 



□ 

Proof of Theorem ICT 

The theorem is immediately imphed by Lemma 1?^ and the following 
result. 

Lemma 4.5. Let (X, $, 7?) be a p-stable sampling model for some p G 
[1, 00] and X = X + A. Let U be the sampling operator for (X, $, 7?) 
and Ua be the sampling operator for (X, $, 7?)- Then \\U — Ua\\ —>■ 
as ||A||oo — > 0. 

Proof. We recall that for any 7 > 0, the function osc^ g on M'^ is defined 
by 

osc^ g{x) = sup \g{x + Ax) — g{x)\. 

|Ax|<7 

From Lemma 8.1 in [4J it follows that if (7 G VFq^, then osc^ g G VF^, 
and II osc^ g\\ ^/i ^ as 7 — > 0. Therefore, by applying Proposition 14.21 
we get 

osc^ $ * 7? G and || osc^ $ * 7?||(vFi)('-x*) ^ as 7 ^ 0, 
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where 



osc^ $ * 7? 



OSC^ (f)^ * fl^ ... OSC^ (j)^ * /i* 



OSC^ (jf * jj} ... OSC^ 0'" * /i* 



For any m e7/ there exist at most {[5 ^\/d\ + 1)'^ samphng points in 
every hypercube [0, l]'^ + m. We set Xm = X n([0, l]'^ + m), m G Z*^, 
and, for each 1 <i <r and 1 < / < t, define the sequence 

6*''(m) := esssup^g[o,i]<* {osc||a||^ (0* * ^Ji!'){x + m)}, m G Z*^. 

Then ||6*''||£i(2d) = || osc||Aj|^ (0'' * w^^ and, hence, 

= l|0SC||A||^ * 7?||(u^i){rxt). 

For 1 < i < r and 1 < / < t we have 



J2 4 ((0^ * - A:)) - (0^ * + 6, - k)) 



where A/" = (^^^v^ + lY'"^- By using Young's inequahty we obtain 

-/V |||C|*0 ll^p(zd) ^ -/V ||C ||^p^^d^||0 11^1 

= Ar*'l|c1|^p(zd)||0SC||Ai|^ 0**/i'f^^i. 

Consequently, 

||f/^''-f/i'l| <Ar||0SC|,A|U 

Hence, 

lit/- t/All < Ar||osc||Aj|^ $*7?ll(vvi)(-*) ^0 as ||A||oo^0, 

and the lemma is proved. □ 

Proof of Theorem [3A0l 

Proof. The proof of Theorem 13.101 is hidden in the proofs of Theorems 
13.11 13.41 and 13.81 In particular, keeping the notation of the proof of 
Theorem 13.11 we have 

ll((/ -9)* 7?)(-^)ll(£P(J))(t) < 2'^A^||C||(<,p(2;d)){r)||$-e||(jyl){r)||7/||(^(Rd))(t). 
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Hence, Theorem 13.101 is true, when ~jj = 'a and X = X + A. Keeping 
the notation of the proof of Theorem 13.4^ we have 

-7?))(X)||(^p(j))(t) < 2'^A/'||C||(£p(2d))(r)||$||(VKi)M||7?-"tt ||(A4(M''))(*)- 
This inequahty imphes Theorem 13.101 when $ = and X = X + 
A. Combining these resuhs with Theorem 13.81 via the standard e/3 
argument we prove the general case. □ 

4.3. Proofs for Section 3.2. 

We begin with an auxihary technical result for the convolution of 
functions with measures. 

Lemma 4.6. Let $ = (0^, . . . , cj)'^)'^ be a vector of continuous functions, 
s>d, and~fx e {Ms(M.'^)Y*\ If < C'q{1 + \x\)-' for all 1 < i < 

T thcTl 

\{<^*lt){x)\<C^{l + \x\)-^■ 
the constants Cq > 0, 1 < i < r, and Ci > are independent of x G W^. 

Proof. For 1 <i <r and 1 < j < t we have 

< I \<p\x-y)W\{y) 



< Q / (l + \x-y\)-^d\fi^iy). 

Since (1 + |m + w|)-' < (1 + |m|)'(1 + for all u,w e M'^, and / > 0, 

we have 



< CU {l + \y\ni + \x\rd\f^^\{y) 
= C^(l + |x|)-^ / (l + |y|W|(2/) 



< ci'\i + \x\r, 

where the last inequality follows from G Ais{^^)- Therefore, 
\i^*lt){x)\<C,{l + \x\)-% 



where Ci = E:=iE;=iC^I''- □ 



Remark 4.1. If {^k}k&d is an s-localized Riesz generator for as 
in Definition 13. 4[ then, by Lemma 14(a) in ^ISj, we have that {^k}keZ'' 
is also an s-localized Riesz generator for Consequently, by 

Lemma [4.61 we have 



(4.16) |($*7i)(x)| <Di(l + |x|)- 

for some Di > independent of x G M."^. 
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Proof of Proposition 13. Ill 

Proof. Let X be a 7?-sampling set for Jt e {M{R'^)Y^\ Then, 

by definition, there exist constants < A2 < B2 < oo such that 

(4.17) A^WfllL^ < < i?2||/||L^,for all / G V'{^). 

Fix Xj G X. Then, for each 1 < i < t, the function gl..: C 
given by glXf) = (/ * is a bounded linear functional on the 

closed subspace V^{^) ofL'^{R^) because \gi^{f)\ < 52||/||l2 for all / e 
V"^($). Consequently, by Riesz representation theorem, there exists 
tpi^ e such that gi^if) = for all / G V^{(!>). It follows 

immediately from (14.171) and Definition 13. II that "^xj = {i'l^, ■ ■ ■ 
is a frame for V^^($). Hence, every / G V^^($) can be recovered via / = 
Eiej(/' ^^•.)^^.' where = {^l^,. . is a dual frame of 

and the series converges unconditionally in \^^($). Since 
(/, ^,^) = (/ * Jt){xj) for all J G J, we get (EJD. □ 

Next, we show that if the generator $ and the measures 7? satisfy 
an appropriate decay condition then the {ji, X)-sampling frame {^z^ } 
obtained above is s-localized. 

Proposition 4.7. Let s > d, ^ e Ws, and 'Jm G (A^,(R'^))W. If X is 
a -sampling setforV^{^), then the (Ji , X) -sampling frame {^z^} is 
s-localized with respect to the Riesz basis {^k}kez<i- 

Proof. Since {$fc}fcgzd is an s-localized Riesz generator for V^($), the 
components of $ satisfy (13.61) . and Lemma [4.61 implies 

vj/y I = |($ * -fjt){x, - k)\ < Ci(i + |x, - k\r, 

for some Ci > independent of j G J and k G Z'^. On the other hand, 
it follows from Remark |4. II that the dual Riesz basis {^fcjfcgz'' is also an 
s-localized Riesz generator for V"^($), and its components also satisfy 
(13.61) . Therefore, using Lemma [4.61 once again, we get 

|($,, vl/^^l = |($ * -it)ix, - k)\ < D,il + \x, - k\r^, 

for some Di > independent of j G J and k G Z'*. Hence, {^z^} 
satisfies all conditions of Definition 13.31 □ 

We conclude this subsection with the proof of the main result of 
section 3.2. 



Proof of Theorem [3A2] 
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Proof. Assume the hypotheses of Theorem I3.12[ By Propositions 13.111 
and \4:.7\ there exists a (7?, X)-samphng frame for \^^($), 

which is s-locahzed with respect to the Riesz basis {$fe}fcez<* and sat- 
isfies 



(/,vl/,..) = (/*7r)(x,), for ^nfev\<i>). 



Moreover, 



/ = * for all / G V'i^). 

Consequently, applying Theorem 10(c) in [15], we get 

/ = Y.^i * 7?)(^.)^.., for all / e v^i^), 

where the series converges unconditionally in V'p{^), 1 < p < oo. 
Moreover, since {^^aj^ jjgj is an s-localized frame with respect to the 
Riesz basis {^fcjfcez'*) then Theorem 10(d) in [15] implies that for each 
1 < j9 < oo there exist < Ap < Bp < oo such that 



ApWfh. < IK/ * 7)(X)||(,.(^))W < BpWfU., for all / G y^(<l>), 
i.e., X is a 7^-sampling set for V^{^) and the theorem is proved. 

4.4. Proofs for section 3.3. 

For the proof of Theorem 13.181 we need the following two lemmas. 
Lemma 4.8. Let the assumptions of Theore'm \3.18\ hold. Then 
\\U*U-UIUa\\ <e{2(3p + e). 



□ 



Proof Since \\U\\ = \\U*\\ and \\U - 



\U* - f/^ll. 



\U*U-UIUa\ 



< 
< 
< 



U*U - U*Ua + U*Ua - U^UaI 
U*{U -Ua) + {U* -UDUaW 



f/*||||f/- f/All - 

U-Ua\\{\\U\\ 

u -uaW mu\ 



\u* -ui\\\\Ua\ 

llf/All) 

■\\U-Ua\\) 



< e{2(3p + e) 



and the lemma is proved. 



□ 



Lemma 4.9. Let the assumptions of Theorem \3.1^ hold. Then < 
V <l, {UIUaY^ exists, and \\{U*U)-^ 



{UIUa)-^ < 



nl{i-v)- 
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Proof. Since {U*U)~^ exists, 

(4.18) uiua = u*u (J + {u*uy^ {U^Ua - U*U)) . 

From (EJD and (ES]) we get that for all C G (£p(Z'^))('-) 

11^*11 (£P(Zd))(r) < \\{U*U) ^C||(£p(2d))(r) < — He'll (^P(Zd))('-)- 

From the above inequahties and Lemma 14.81 we have 

\\{U*U)-^{UIUa-U*U) II < ||(f/*f/)-i||||f/^f/A - f/*f/|| 
< (2/3, + e) < ^(-/3, + V^T^) (2/3, - /3, + V'^T^) = 1. 

Hence, = \e (2/3„ + e) G (0, 1). To simphfy the notation, we define 

M ■= U*U, Ma := U^Ua, and := (f/*f/)-i {U*aUa - U*U) . 

Since ||A^|| < u < 1, then (/ + A^)~^ exists and is given by the Neumann 
series 

oo 

(J + A^)-i = ^(-1)'?A^9. 

q=0 

From ( ]4.18p we obtain 

(4.19) M^^ = [M{I + N)y^ = (/ + N)-^M'\ 

Therefore, M^^ = {UlUAy^ exists. 

Now we need to give an upper boundfor ||A^-i-M^^||. Using (gli) 
we obtain 

M-^ - M^^ = N{I + N)-^M-\ 

Consequently, 

\\M-^-MX^\\ < ||iV||||(/ + A^)-i||M-^|| 

ll^^ll „ , r 1 „ ^ 9 

- 1 _ ||Ar|| " II - 1 _ I, 'p ' 

(4.20) 

and the lemma is proved. □ 
Proof of theorem 13.181 
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Proof. Using the notations from Lemmas 14.81 14.91 and the previous 
proofs, we get 

\\{u*u)-^u* - {uiUAy^uiw = \\M-'^U* - M^^UIW 

= \\M-^U* - M-^Ul + M-^Ul - M^^UIW 
= \\M-\U* - Ul) + (M-i - M^^)Ul\\ 
< ||M-^||||?7* - UIW + \\M-^ - M^^WWUIW 

Vp\ l-v J 

□ 

Proof of Theorem [3A9l 

Proof. Let f/^ be the samphng operator for a perturbed samphng model 
{X + A,e,^). Let also C e / = EkeZ'^ Cl<l>k, and g = 

Ek&'i ClQk- Then 

\\R{g * + A) - /ll^, < M, \\{{WU)-'WUaC - C) ||^, . 

It remains to apply Theorem 13.181 to finish the proof. □ 

Proof of Theorem [QUI 

Proof. Assume the hypotheses of Theorem 13.201 From Theorem 13.121 
we know that, in this case, the sampling model (X, 7?) is p-stable 
for every p G [1, oo]. Hence, in view of Theorem I3.19[ the only thing 
that we need to prove is that the operator U*U is invertible for all 
p e [1, oo] and not just for p = 2. 

Taking into account that for each 1 < i < r and I < I <t the entries 
of the matrix of the operator f/*'' satisfy 

\iU^\k\ = 1(0* * f^')ix, -k)\< Cil + \xj - A;|)-^ 

for some Ci > independent of j G J and A; G Z"', it follows from 
Lemma 3 in [15j that the matrix of U defines a bounded linear operator 
from ^ (^P(J))(*) for all 1 < p < oo. Hence, U* is also well 

defined as a bounded linear operator from {i^^J))^^^ and, 
therefore, U*U : (£p(Z'^))(^) ^ (£p(Z<^))('') is a well defined and bounded 
operator for all 1 < p < oo. On the other hand, since the operator 
U*U is invertible on (£^(Z'^))('') and its components (M*'')j^fc, 1 < i < r, 
I < I < r, satisfy a decay condition 

m'''),,k\<c2{i + \x,-k\r% 

for some C2 > independent of j G J and k G Z'', then Jaffard's 
Lemma (see Theorem 5 in [15]) implies that {U*U)-^ : (£2(Z'^))(^) 
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(i'^ (Z'^))^^^ is also a bounded linear operator defined by a matrix sat- 
isfying the same off-diagonal decay condition as U*U. Consequently, 
using Lemma 3 in [15] once again, we get that the matrix of {U*U)~^ 
defines a bounded linear operator on (£p(Z'^))^'''^ for all 1 < p < oo. 
The theorem is proved. □ 
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